Kramers-Kronig constrained variational analysis of optical spectra 
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A universal method of extraction of the complex dielectric function e(ui) = £i(w) + 162(01) from 
experimentally accessible optical quantities is developed. The central idea is that 62(01) is param- 
eterized independently at each node of a properly chosen anchor frequency mesh, while ei(oi) is 
dynamically coupled to 62(01) by the Kramers-Kronig (KK) transformation. This approach can 
be regarded as a limiting case of the multi-oscillator fitting of spectra, when the number of oscil- 
lators is of the order of the number of experimental points. In the case of the normal-incidence 
reflectivity from a semi-infinite isotropic sample the new method gives essentially the same result 
as the conventional KK transformation of reflectivity. In contrast to the conventional approaches, 
the proposed technique is applicable, without readaptation, to virtually all types of linear-response 
optical measurements, or arbitrary combinations of measurements, such as reflectivity, transmission, 
ellipsometry etc., done on different types of samples, including thin films and anisotropic crystals. 



I. INTRODUCTION 

The linear macroscopic optical characteristics of mate- 
rials, such as reflectivity, absorption, penetration depth 
and others, are largely determined by the complex dielec- 
tric function e(w) = 61(01) +{62(01)^. Typically, one has to 
extract 61(01) and 62(01) of the studied sample by the in- 
version of measured optical spectra. These quantities are 
not independent, but comply the Kramers-Kronig (KK) 
relationa 2 i 3 i 4 i 5 : 
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■ where <jdc is the DC conductivity, and the symbol p 
denotes the principal value integral. 

When talking about experimental ways to extract of ei 
' and 62 , one should distinguish between the ones that rely 

on the KK relations and the ones that allow the determi- 
' nation of both quantities directly from the experiment. 

A well-known example of the first type of techniques is 

■ the KK analysis of the normal-incidence reflectivity R(oi) 
spectra 6 -. The idea of this method is the following. Equa- 

■ tions QJ and @ cannot be directly applied, since both 
6\ and 62 depend on the unknown phase 9 of the complex 
reflectivity 



1 + V~6 



i?exp (id). 
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A very helpful 'trick' is to write down a relation similar 
to (J5J) with respect to the complex function lnr(aj) = 
lny/R(u) +i6(oi): 



e(ui) = --P f 
n Jo 



In y/R(x) 



dx + 9(0) 



(4) 



and use it to compute 6(01) from R(oi). The reflectiv- 
ity has to be extrapolated outside the experimentally 
accessible frequency range, which may give rise to a 
large uncertainty, if the experimental range is not suf- 
ficiently broadA. Having found the reflectivity phase, 
one can restore the dielectric function by inverting (J3J: 
6 = (1 — r) 2 /(l + r) 2 . Thus, the KK relation is an ad- 
ditional condition, which 'miraculously' helps to get two 
spectra, 61(01) and 62(01), out of one, R(o>). 

Unfortunately, the approach based on equations like 
J3J, is not universal 8 -. The KK relations are derived for 
the so-called response functions^, which have regular an- 
alytical behavior in the upper complex frequency half- 
plane due to the causality principle. In the case of the 
function lnr(w), a necessary condition is that r(oi) =/= 
for Im(o;) > 0. One can show, that it holds for the nor- 
mal angle of incidence, but breaks down, for instance, for 
off-normal reflectivity, if the angle of incidence is large 
enough^. The same problem appears, when the sample 
is a thin film on a substrate and in many other cases. 
A number of analytical and numerical algorithms were 
put forward to extend the KK method to off-normal 
reflection 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 , transmission 16 ! 17 ! 18 ! 19 , attenu- 
ated total reflection (ATR) spectra 13 ' 20 ! 21 ! 22 as well as to 
optical data on layered samples 1 ^* 1 ^* 2 ^ 4 * 2 ^ 2 ^ 7 . and low- 
symmetry crystals^iS 9 . and other 'non-standard' situa- 
tions. Although the proposed techniques, taken together, 
certainly cover a substantial scope of applications, it is 
desirable, from the practical point of view, to have a sin- 
gle universal approach, which can be used in most of the 
currently used and future optical experiments without 
major readaptation. 

The least-square fitting of spectra using simple model 
dielectric functions, such as Drude-Lorentz oscillators, 
also called the dispersion analysis 3 ^, is practically more 
robust^ 6 ^ than the mentioned above KK-based tech- 
niques. However, the range of considered dielectric func- 
tions is artificially restricted to a certain model, which 
makes it hard to recover, without distortions, all signif- 
icant spectral details from the experimental data. On 
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the other hand, the introduction of new parameters to a 
model (for example, the addition of extra Lorentz oscil- 
lators), makes the functional form for the dielectric func- 
tion more 'flexible' and, in a sense, less model-dependent. 

The method of the KK- constrained variational fitting, 
developed in this paper, is a 'hybrid' of the two ap- 
proaches. It inherits the robustness of the least-square 
fitting and the ability of the KK analysis to extract the 
full spectral information. In simple terms, it is a limiting 
case of the multi-oscillator modelling, when the number 
of oscillators is so large (of the order of the number of ex- 
perimental points), that the analysis becomes essentially 
model-independent. Since the KK relations (11 1211 are au- 
tomatically preserved in the used functional form of the 
dielectric function, the proposed approach can also be 
regarded, as an extension of the standard KK methods. 

The basic idea of the KK-constrained variational fit- 
ting was already introduced in Ref^S for the special 
case of the normal-incidence reflectivity of low-symmetry 
(monoclinic and triclinic) crystals, where the standard 
KK method fails. The KK relations for the diagonal and 
off-diagonal components of the reflectivity tensor were 
used, which limited the applicability of the technique. 
Here we overcome this limitation by referring directly to 
the KK relations for the dielectric function. 

The rest of the manuscript is organized as follows. In 
Section[n]we outline the general problem targeted by the 
method. Section Hill discusses the variational fitting as a 
limiting case of the usual data fitting with a large number 
of oscillators. In Section Hvl we describe the technique in 
details. Section Ivl contains a number of application ex- 
amples. Finally, in Section IVII we summarize the results 
and overview the possible future developments. 

II. PROBLEM FORMULATION 

We begin with a formal problem definition. Sup- 
pose, a set of experimental datapoints {co^ xp , q^ xp , Sq^ xp } 
(i=l,...,N exp ) is given, where u>^ xp is frequency, q exp is an 
optical quantity of a certain type (reflectivity, transmis- 
sion, penetration depth, an ellipsometric angle etc.) and 
5q^ xp is an experimental error bar of q exp . We assume, 
that all optical quantities can be expressed via real and 
imaginary parts of the material dielectric function e(uj): 

c' c = /^i(^r)^r p ))- (5) 

Different datapoints can be of different types, therefore 
we keep the index i in fi, Using the transformation 
(which can be symbolically written as {ei} — 1 = 
KK({e 2 })), we get: 

C' C (fe}) = M1 + KK{e 2 }(u*">)Mu? P ))- (6) 

Note that qf alc depends not only on e 2 at ijf xp ', but also 
on £2(0;) at all other frequencies. 

It is natural to look for a dielectric junction e{uj), which 
satisfies the KK relations and provides the best fit of all 



experimental datapoints: 

X ({£2}) = 2_, J^p j -'•nun. (7) 

This relation formalizes the least-square approach, which 
is justified in the case of the normal error distribution 32 . 

The minimization problem (JJJ is not fully defined un- 
less we specify the functional space for e 2 (u>). In princi- 
ple, one may consider all physically acceptable functions 
e 2 (u>) (0 < lo < 00), which are non-negative and inte- 
grable (i.e. compatible with the f-sum rule2£): 

POO 

i xe 2 (x)dx < 00. (8) 
Jo 

However, in such a broad functional space the minimiza- 
tion problem becomes ill-posed^ k In particular, multiple 
solutions for e 2 {oj) may exist, which look very different, 
but minimize (J7J equally well. However, if the set of data- 
points is sufficiently broad and dense, only one (or maybe 
few) solutions will show a regular behavior, while the ma- 
jority of them will be marked by 'suspicious' features, for 
instance, sharp spikes, discontinuities or oscillations with 
a period smaller than the distance between experimental 
frequencies. 

It suggests that a proper way to regularize this ill- 
posed problem is to restrict somehow the functional space 
e 2 [uS). It can be done using some parametrization, for 
example, the decomposition of e 2 (u>) in a finite basis of 
fixed reference functions. An optimal parametrization 
should exclude the mentioned irregular features, while 
being flexible enough to fit all important spectral details. 
Without rigorous mathematical proof, but rather based 
on practical tests, we argue that the optimal number of 
parameters should be of the order of but not larger than, 

One should note that it is more common to use model 
dielectric functions with a small number of parameters, 
much less than N exp . In this case, the purpose of the 
analysis is to get a physical insight using the measured 
spectra directly. For instance, the validity of a certain hy- 
pothesis about the studied compound can be tested and 
the values of some meaningful parameters (such as the 
plasma frequency, optical transition energies etc.) can be 
found. On the contrary, the purpose of parametrization 
in the model-independent method, proposed in this paper, 
is not to get the values of parameters (which may have 
almost no physical meaning, taken alone), but rather to 
extract the complex dielectric function by the inversion 
of the available data. Thus, it is not surprising that the 
number of parameters becomes very large. The mini- 
mization of Q in this case is essentially a variational 
problem. The obtained in this way dielectric function 
can be used later on for a more specific examination. 
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III. DRUDE-LORENTZ FITTING AS A 
PARAMETRIZATION 



An important example of parametrization of dielectric 
functions is the Drude-Lorentz (DL) oscillators formula 



N 



eiuj 



k 



1 U lk 



(9) 



Each oscillator is represented by a Lorentzian with three 
adjustable parameters: the oscillator frequency took, the 
linewidth jk and the plasma frequency LU P} k- The un- 
bound oscillators (conductivity electrons) are described 
by Drude terms, where wo.fc = 0. The parameter is 
the contribution from the higher- frequency oscillators. 

Usually the number of oscillators is very limited, com- 
pared to N exp . As a result, the obtained spectra of e(to) 
are 'unrealistically' smooth with some subtle, but po- 
tentially important, spectral information being lost. Ac- 
cording to previous arguments, on should take a very 
large number of oscillators (of the order of N exp ) in or- 
der to extract the information model-independently. It 
was argued (see Refsi2§i2i) that the fitting of the most 
important details of spectra (other than statistical noise) 
is not much different from the KK transformation. The 
physical meaning of individual oscillator terms may not 
be always clear, however, the meaningful result is the 
complex value of e(w) itself. 

In practice, the oscillators are usually introduced one 
by one and all Lorentzian parameters are allowed to 
change. Although this strategy works very well for small 
N (~ 10 - 20), it becomes cumbersome and badly con- 
trolled as the number of oscillators grows further. In 
particular, as one adds more and more Lorentzians to 
the model function JJjJ, it becomes less and less clear 
how to guess the initial parameter values. If all parame- 
ters are free to change, it often happens that the fitting 
process becomes numerically unstable. To avoid it, one 
has to deliberately fix some parameters. It makes the fit- 
ting procedure rather tricky and ambiguous, preventing 
the routinely operation. In the next section we elaborate 
another fitting scheme which is more attractive from the 
practical point of view. 



IV. METHOD 

Now we proceed with the detailed description of the 
KK-constrained variational fitting algorithm. Although 
it can be concisely described in a hand-book like man- 
ner, we prefer to show its 'evolution' starting from the 
described above DL fitting scheme. In this way, one can 
easily see that the new technique is, in principle, a lim- 
iting case of the multi-oscillator fitting. 



A. Setting up the oscillators 

Let us consider a dense enough mesh of anchor fre- 
quency points u>i,..., ujn- It must not the same as the set 
of experimental frequencies Lo e ^ p , introduced in 

Section |nl although they may coincide. We shall discuss 
below the optimal choice of the anchor mesh. 

It is possible to set the number of oscillators in ©, 
their frequencies and linewidths in a rather logical and 
simple way. To be 'flexible enough', the fitting function 
J5J should be able to deliver spectral weight everywhere in 
the spectral range It can be done by putting one separate 
oscillator at every anchor frequency: ujq^ — uji . Further- 
more, we would like the i-th oscillator to be narrow and 
provide spectral weight only in the closest neighborhood 
of u) i . It suggests to set its width to a fixed value of the 
order of the distance between adjacent frequency points, 
for instance ji = (to>i+i — cJi_i)/2. Figure^a) shows the 
real and the imaginary part of the i-th Lorentzian: 



LO 2 ) 2 + jfoJ 2 
,2 



4fH 



(10) 



(11) 



(uj 2 -UJ 2 ) 2 +-/ 2 uj 2 

In this way we get a dense 'forest' of narrow 
Lorentzians, sitting at fixed anchor frequencies. Only 
the oscillator strengths Si are left adjustable in order to 
set up the profile of the 'forest'. 
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FIG. 1: Dielectric functions with a Lorentz (a), a triangular 
(b) and a rectangular (c) lineshapes, centered at uJi and having 
the linewidth 74 = (wj+j — Ui_i)/2. 
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B. Optimizing the oscillator lineshape 

Because of extended low- and high-frequency tails, 
which interfere with other oscillators, the Lorentz line- 
shape (Figure ^a)) is not very suitable to represent the 
'local' spectral weight near The problem becomes 

evident from the following example. Imagine, that there 
is a low- frequency gap uj g in the spectrum: t2{uj) = for 
uj < ujg. However, each Lorentzian that we have to put 
above uj g is going to contribute a non-zero value of €2 
below ujg because of its low-frequency tail. This makes 
a good fitting quality with Lorentzian oscillators hardly 
achievable. 

Clearly, it is better to have a lineshape function, which 
is non-zero only inside a small region, adjacent to Wj. A 
good candidate is a triangular profile, shown on Figure 



A 
£ 2,i 



(U - U)i-i)/(u>i - Wi-l) , CJj_i < UJ < UJi 

(uj) = <( (oj i+1 - uj)/ (u>i+i - uji) , UJi < uj < uj i+1 . 

, otherwise 

(12) 

The real part can be obtained by the KK transformation 



A / \ 1 
7T 



g(w,Wi-i) (uj i+ i - uJi-x)g{uj,uJi) 



(UJi 



-l)(uJi+l - UJi) 

g{uj,uj i+ i)' 



H+i 



(13) 



with g(x, y) = (x + y) In \x + y\ + (x - y) In \x - y\. 

Exactly as we did it before, the variational dielectric 
function can be constructed as a linear superposition of 
triangular functions sitting at fixed anchor frequencies 



OJ 

to" 




Frequency 



FIG. 2: An example of the variational dielectric function 
€ va r(uj), described in the text. The imaginary part of 
superposition of several triangular functions e^ i (u>), centered 
at anchor frequencies. ei,„ar(k>) is the exact KK transform of 

£2,var(uj). 



of which quadrature formula to use for the numerical in- 
tegration in the KK transformation. The triangular line- 
shape shape corresponds to the trapezoidal integration. 
One can also use a rectangular shape of e^^uj) (Figure 

^c)). However, it is inferior to e^(w), since the discon- 
tinuity of e^i^uj) would result in artificial divergences of 
e^i^uj) at the rectangle edges. 

C. Discrete Kramers-Kronig transform 

Since we deal with the values of the dielectric func- 
tion at a fixed mesh, it is convenient to represent its 
real and imaginary parts as vectors: eu = e\{uii) and 
e 2 ,i = £2(^1)- Then, using formulas lfhf|) and (fT3"|) . the 
KK transform Q can be represented in a matrix form: 



N-l 



(14) 



TV 



(15) 



The coefficients Ai are tread as free parameters. We ex- 
cluded i = l and i = N , from the summation by defining 
Ax = A N = 0. 

The variational function is shown schematically on Fig- 
ure[21 The imaginary part of e var {uj) is a piecewise curve 
passing through points {uji, Ai} (i — 1 . . . N). Thus, the 
adjustable parameters have very straightforward mean- 
ing: they are the values of €2 at anchor frequencies Wi. 
Obviously, any shape of C2{uj) of a real material can be 
reproduced with a proper set of parameters Ai, provided 
that the frequency mesh is dense enough. 

At first glance, a triangular lineshape is rather 'patho- 
logical' for an oscillator. However, we are going to ascribe 
the physical meaning not to the individual oscillators, but 
to the total dielectric function (H) , composed of many 
of them. It makes us free to choose a lineshape, which 
is the most convenient from the computational point of 
view. Even though this choice is rather arbitrary, it is 
not crucial for the method. It is equivalent to a decision 



where 



K, 



A 1 \ 



(16) 



Keeping in mind the analogy with the discrete Fourier 
transform, we can call (|15JI a discrete Kramers-Kronig 
transform (DKKT). The matrix Kij is a discrete analog 



of the KK kernel K(x,uj) = (2/n)x{x 2 



The exact expression for the DKKT depends, of course, 
on the chosen oscillator lineshape. However, for > 2 
all lineshapes give almost the same value of K^. 



D. The problem of extrapolations 

A serious shortcoming of the variational function l|14|) , 
taken alone, is that it does not take into account any spec- 
tral weight beyond the frequency range uj\ < uj < wjy. 
Due to the non-local character of the KK relation (Q, 



5 



the non-zero wings of £2 at higher and lower frequencies 
are essential to calculate £1 inside this region. Therefore, 
some assumptions about the high- and low-frequency be- 
havior of £2 should be taken. 

The same problem of extrapolations exists in the usual 
KK analysis of reflectivity, which does not have any uni- 
versal solution. Nevertheless, an 'educated guess' can 
often be done, based on some reasonable physical model. 
We argue that the following 'two-step' fitting scheme can 
address the issue of the extrapolations. 

At first, the experimental spectra can be fitted in a con- 
ventional way by some formula-defined dielectric func- 
tion e mo d(w) with a limited number of parameters. It 
can be a DL function @ or a more appropriate physical 
model. If e mo d(uj) describes the major features of ex- 
perimental spectra then we can assume that it gives an 
approximately correct frequency dependence outside the 
considered spectral range, even though some fine spec- 
tral details may not be fitted very well. Of course, in 
order to have a better confidence in the extrapolations, 
the experimental spectral range has to be expanded as 
much as possible, which is also true for the conventional 
KK analysis of reflectivity. 

The next step is to add the variational function Ijl4(l : 

etotal(u) = e mod (u) + e var {uj), (17) 

to 'fix' the parameters of £ mo d(w), and to perform the 
final data fitting with only parameters Ai of the varia- 
tional function (|14|) kept adjustable. The low- and the 
high-frequency spectral weights are now accounted for 
by e mo d{uj). In this case the 'flexible' function e var {u>) 
acts as a small correction, intended to fit the fine spectra 
details, unaccessible for the 'stiff' e mo d{oj). Therefore, 
the initial guess values Ai can be set to zero. Note that 
even though e2,totai{^) ought to be positive, the sign of 
t2,var{u) is now arbitrary. 

E. Computational details 

Once the set of datapoints, the physical models and 
the anchor mesh are specified, one can minimize the chi- 
square functional Q using of the standard non-linear 
minimization algorithms. In particular, the Levenberg- 
Marquardt (LM) method^ has proven to be very effi- 
cient, when the number of parameters is large. An im- 
portant advantage of the LM algorithm is that it involves 
the calculation of derivatives of the model spectra with 
respect to parameters using explicit analytical formulas. 

The optimization of the mesh of anchor frequencies 
lui,..,lun is rather important for the correct work of the 
variational fitting routine. On one hand, the mesh should 
be dense enough to enable the variational function (|14|) 
to fit all important spectral details. On another hand, N 
should not be too large. The first reason is that the cal- 
culation time grows quickly as a function of the number 
of parameters to be adjusted. The second reason, dis- 
cussed in Section m is that an excessive 'flexibility' may 



cause irregular behavior of the resulting functions, such 
as fast spurious oscillations. According to our experience, 
for an optimal mesh, N is two times smaller, than N exp , 
while every anchor point corresponds to at least one ex- 
perimental point, contributing to formula J7J). On an 
everyday-use computer, one can handle up to few thou- 
sands parameters within reasonable time bounds, which 
is sufficient in most cases. 

Another way to avoid physically meaningless numerical 
oscillations would be to add extra 'regularization' terms 
to the total chi-square functional JTJl, intended to make 
unfavorable those solutions, which have short-period os- 
cillations. In this paper we are not going to discuss the 
details of possible regularization schemes. This work is 
in progress. 

If the mesh optimization and other precautions fail to 
get rid of a numerical instability, it may signal that the 
systematic error bars are too large, or an inappropriate 
physical model is used. Such 'inconsistency alert', which 
is missing in the usual KK method, is another advantage 
of the proposed technique. 



V. APPLICATION EXAMPLES 

Now we turn to examples, showing the applications 
of the method. All calculations were performed using 
Rcffit software^, designed for the fitting of optical spec- 
tra with various physical models as well as the KK- 
constrained variational dielectric functions. As an exper- 
imental input, we shall use either really measured curves, 
or computer-generated test spectra. 



A. Normal-incidence reflectivity 

Let us first apply the variational fitting method to a 
normal-incidence reflectivity spectrum and compare it 
with the results of the usual KK analysis 6 . Figure |3fa) 
shows the far-infrared room-temperature reflectivity of a 
high-T c superconductor Bi 2 (Sr,Y) 2 CaCu208 (shortly Bi- 
2212) for the electric field along the c-axis (perpendicular 
to the planes Cu02). The spectrum is dominated by op- 
tical phonons, since the interplane electronic transport is 
almost blocked. 

We applied the usual KK analysis, extrapolating the 
reflectivity by a constant value at high and low frequen- 
cies, because the conductivity of Bi2212 along the c-axis 
is very small. 

In order to perform the variational fitting, we took 
a mesh of 700 anchor frequencies, logarithmically dis- 
tributed between lj\ = 60 cm -1 and u>n = 1500 cm . 
The resulting real and imaginary parts of e(uj) are shown 
in Figure [3fb) and (c). One can see that the two tech- 
niques give almost identical results in this reference case, 
when the conventional KK transform of reflectivity is 
known to be operational. 



6 



0.8 




Wavenumber (cm' 1 ) 

FIG. 3: The KK-constrained variational fitting vs. the usual 
KK analysis, applied to the normal-incidence reflectivity of 
Bi-2212 for E || c, at room temperature (a). The two methods 
give almost indistinguishable £i(o;) (b) and £2(w) (c). The 
difference between the results, multiplied by 10, is shown for 
comparison. 

B. Grazing-incidence reflectivity 

The measuring of reflectivity R p (uj) at a grazing angle 
of incidence (close to 90°) for the polarization of light in 
the plane of incidence (p-polarization) is often preferable 
for highly reflecting metallic samples, since the measured 
absorption A = 1 — R p is more sensitive to the sample 
conductivitjiSi, than it is for the normal reflectivity. We 
are going to show that the variational fitting method is 
workable in the case of the off-normal incidence reflectiv- 
ity. Here we use an artificial example, where the 'true' 
dielectric function is known beforehand. 

The 'true' dielectric function e(uj) = £i(w) + 
47T«(Ti(w)/w, which is shown in Figure 0] (b) and (c), 
was simulated by a DL model 10 with 10 randomly put 
Lorentz oscillators, one Drude peak and eoo=6. The 'ex- 
perimental' reflectivity for the p-polarized light for the 
angle of incidence 8=80° was calculated by the Fresnel 
formula-: 



e cos 6 — \/ e — sin 2 6 

p = / , | • 

e cos 9 + v e — sin o 

Some statistical noise of 1% was added (red solid line on 
Figure0](a)) in order to mimic a real-life situation. We 
verified that the usual KK analysis does not work in this 
case: the formula (@J, when applied to R p , gives a totally 
wrong reflectivity phase spectrum. 



The obtained in this way R p (lu) was taken as an input 
for the KK-constrained variational fitting procedure. At 
the first stage the rough fitting was performed with only 
3 oscillators included to e mo d{uj) (green dashed line on 
Figure 0(a)). After that a variational function e var (uj) 
was added to e mo d{u) and the fine fitting of R p (uj) was 
done. We took a mesh of 500 anchor frequency points, 
logarithmically distributed between 50 cm -1 and 10000 
cm -1 . The resulting ei(ui) and <Ji(u>) (Figures 0] (b) and 
(c), red line) are very close (apart from the noise) to the 
'true' ones (Figure 0(b) and (c), blue line). Of course, 
such a good agreement is only possible, when the fre- 
quency range of the input reflectivity spectrum is large 
enough. 




Wavenumber (cm" 1 ) 

FIG. 4: The KK-constrained variational fitting of a grazing- 
incidence reflectivity spectrum for 0=80° (a, red solid curve) , 
simulated, using 10 Lorentz and 1 Drude oscillators as de- 
scribed in the text. The rough 3-oscillator fit is shown by the 
dashed green curve. The resulting £i(w) (b) and (Ji{uj) (c) are 
very close to the 'true' spectra. The difference between the 
results, multiplied by 10, is shown for comparison. 

An example of the application of this procedure to a 
really measured grazing-incidence reflectivity spectrum 
(in combination with ellipsometry at higher frequencies) 
can be found in Ref—i. 



C. Reflectivity + transmission 

Now we consider an example, when two measured spec- 
tra of different types are analyzed simultaneously. Fig- 
ure Ela) shows reflectivity and transmission spectra of 
Lai.88Sro.i2Cu04 (LSCO) in the far- and mid-infrared 
range at room temperature^. The electric field was al- 
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ways parallel to the c-axis. The near-normal incidence re- 
flectivity R(u>) was measured on a thick non-transparent 
sample. The transmission spectrum T(u>) was taken on a 
thin slab of the same material with a thickness d—23 fim. 
Both spectra can be expressed in terms of the complex 
dielectric function e(w)&: 



R 



\r\ 2 ,T = 



(1 - r 2 )t 



1 - r 2 t 2 



(18) 



where 



t = exp 



The formula for T takes into account multiple internal 
reflections, which give rise to Fabry-Perot interference 
fringes. Unfortunately, the fringes are partially sup- 
pressed in the experiment, because the sample was not 
exactly flat-parallel. This was taken into account by the 
averaging of T for a certain distribution (about 10%) of 
the sample thicknesses. 

Figure EJb) shows the optical conductivity <7i(a;), ob- 
tained by two different techniques: the standard KK 
transformation of reflectivity (not considering the trans- 
mission spectrum) and the KK constrained variational si- 
multaneous fitting of reflectivity and transmission (which 
we call RT+KK method 38 ). The c-axis conductivity is 
featured by a slowly varying small electronic conductiv- 
ity, determined by an incoherent electron transport along 
the c-axis, and the two strong phonon peaks at 235 and 
495 cm -1 . 

One can notice, that the two techniques give similar 
results in the phonon range, but significantly deviate at 
lower and especially at higher frequencies (watch the log 
scale!), where the conductivity is small. It is not surpris- 
ing, since the reflectivity in this case depends mostly on 
the ei and is not very sensitive to a±. Thus, the error 
bars of the latter are very large. On the other hand, the 
transmission is directly related to the absorption, and, 
correspondingly, to the conductivity. Therefore, it is re- 
ally worthwhile to use both R(uj) and T(u>) in the anal- 
ysis. 

The advantage of using the KK constraint here is not 
immediately obvious, since the two values R and T can 
be in principle converted to t\ and o\ at every given fre- 
quency (RT-method). However, this conversion works 
well only if the value of transmission is not too small. 
In our case the transmission vanishes between 200 and 
600 cm" 1 (in the phonon range), where we have to re- 
sort to the KK relations. Using two different techniques 
in different ranges would not be very practical. Alter- 
natively, the RT+KK method provides a seamless con- 
ductivity spectrum with small error bars in the whole 
spectral range. In addition, the using of the transmis- 
sion spectrum outside the phonon range in the RT+KK 
technique suppresses the uncertainties of the KK analy- 
sis due to the high- and low- frequency extrapolations. A 
similar idea of 'phase anchoring' is discussed in the next 
example. 
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FIG. 5: Reflectivity (bulk sample) and transmission (thin 23 
sample) of LSCO for E||c at room temperature, (b) The 
conductivity spectrum obtained by the RT+KK analysis, as 
described in the text. 



D. Low- frequency reflectivity + high-frequency ei 

and £2 

It is known, that the accuracy of the KK analysis can 
be greatly improved if both real and imaginary parts of 
e(u>) are measured independently in selected frequency 
points 3 ^^ or in a certain interval^i^. Here we consider 
a situation, when one can obtain ei(w) and ei{uS) at high 
frequencies, but only the reflectivity R(uS) can be mea- 
sured at low frequencies. 

In order to extract the conductivity at low frequencies, 
we have to use the KK relations. It is desirable to take 
as much profit as possible from the knowledge of both ej 
and 62 at high frequencies. Therefore, it would be rather 
wasteful, for example, to calculate R from t\ and £2 at 
high frequencies and use it only to extend the range for 
the usual KK transformation On the other hand, the 
simultaneous KK-constrained fitting of the low-frequency 
reflectivity and the high-frequency e\{u)) and €2 (c^) ? 1 '?7i 4 ? 
allows one to use all available information, in particular, 
to 'anchor' the phase of the reflectivity. 

An example is shown on Figure The data were 
taken on a polycrystalline sample of MgB2, which is a 
good metal. Even though the magnesium diboride is an 
anisotropic compound, in the spectra analysis it was as- 
sumed that the powder sample behaves as an optically 
isotropic one, described by an effective dielectric function 
e{ui). The reflectivity at a nearly normal incidence was 
measured in the range 20 - 6000 cm" 1 (Figure [(J a)), and 
the ellipsometry at an angle of incidence 0=80° provided 
ei(w) and e 2 (J) at 6000-37000 cm" 1 (Figure fflb)). 

The variational fitting was done in the same way, as 
described before. The anchor mesh covered both spec- 
tral ranges. The quality of the variational fit is very 
good, giving confidence in the resulting optical conduc- 
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FIG. 6: The KK-constrained variational fitting of the re- 
flectivity in the far- and mid-IR (a) and and €2(ui) in 
the near-infrared to UV ranges of a polycrystalline sample of 
MgB2 at room temperature. The panels (c) and (d) show the 
obtained effective optical conductivity. 



tivity spectrum (Figure Efc-d)). 



E. Reflectivity of a thin film on a substrate 

The next example shows that the method can also be 
used to analyze optical data of layered samples. We shall 
consider a simple but practically very important case 
of an optically thin film grown on a substrate. If the 
penetration depth is comparable, or larger than the film 
thickness then the substrate contributes to the overall 
optical response. One has to take it into account, when 
extracting optical properties of the film from the reflec- 
tivity spectra™. 

We generate dielectric functions of the film e/(w) and 
the substrate e s (u>) in the infrared range using certain 
models in order to check later on the accuracy of the final 
result. For the substrate we took a model that gives op- 
tical properties close to the ones of SrTiOg^i, often used 
as a substrate material. The reflectivity R(uj) and optical 
conductivity <J\{ui) of the substrate are shown by green 
curves in Fig. [7| (a) and (b) respectively. The infrared 
spectrum of SrTiOs is dominated by strong phonon lines 
at 100, 180 and 550 cm -1 . To simulate the properties 
of the film, a model was taken, which gives a metallic 
conductivity with a plasma edge at around 10000 cm -1 
and some additional double-peak structure at about 300 
- 400 cm -1 . The corresponding R{lo) (for a bulk sample) 
and ai(to) are presented in Fig. [7|(blue curves). 

The normal-incidence 'experimental' reflectivity of the 
film+substrate system was calculated using formula™ 



r f + t 2 r fs 



l + t 2 i 



f r fs 
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FIG. 7: A test example of the KK-constrained variational 
fitting of a reflectivity spectrum measured on a thin film (200 
nm) grown on a substrate with known optical properties. The 
test dielectric functions of the film and the substrate were 
generated as described in the text. The panel (a) shows the 
generated reflectivity of a thick film (blue), the one of the thin 
film (red), used as an input, with a noise of 0.3% added, and 
the one of a bare substrate (green). The black curve is the 
initial rough fit of the 'measured' reflectivity. The panel (b) 
indicates the exact optical conductivity of the film (blue), the 
extracted one by the method (red) and the conductivity of 
the substrate itself (green). On both panels the arrows mark 
the optical peaks, which come from the film. 



where 



r S = 



1 



r fs = 



^1+ 



exp yi — y/ej d 



(19) 



for the film thickness d of 200 nm (the red curve in 
Fig[7fa)). Some artificial noise of 0.3% was added, in or- 
der to check the influence of the statistical errors on the 
result. One can see that the reflectivity spectrum of the 
thin films is very different form the one of a thick sample. 
The structures in the far-infrared range are caused not 
only by peaks in the film conductivity (shown by arrows), 
but also by the phonon lines of the substrate. 

As usual, the fitting process was done in two steps. Ini- 
tially, the 'experimental' spectrum was fitted with a very 
limited model e mot j(w). The corresponding reflectivity 
is shown by the black curve in Fig[7{a). Secondly, the 
variational fitting was performed, using e var , defined on 
about 500 anchor frequency nodes (one-half of the total 
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number of experimental points). 

The resulting optical conductivity is shown by the red 
curve in Fig[7{b). The agreement with the exact conduc- 
tivity is rather good. The double-peak structure, which 
was not described by the initial limited model, is repro- 
duced quite well. One can see that the experimental noise 
has little effect at high frequencies, but is strongly am- 
plified below 200 cm . This is hardly surprising, since 
the KK method is known to produce large error bars for 
metallic samples at low frequencies, where the reflectivity 
is close to f . 



VI. SUMMARY 

We introduced a new technique of the extraction of the 
real and the imaginary parts of the dielectric function 
from optical spectra. The idea is that all available data- 
points are fitted simultaneously, using a dielectric func- 
tion, which has a fully flexible (variational) imaginary 
part C2(w) and the real part ei(uj), dynamically coupled 
to the 62 (ui) via the KK relation. In simple cases the 
method works as well, as the usual KK transformation 
of reflectivity 

It was shown that the scope of applications of the vari- 
ational KK-constrained fitting is extremely broad. It ex- 
tends to various optical techniques, or combinations of 
them, as well as to different types of samples. The com- 
putational time required for a typical KK-constrained fit- 
ting session is somewhat longer, compared to the usual 
KK transformation, but it is yet well within reasonable 
limits. It is important, that the method needs almost no 
readaptation for different experimental situations, saving 
a lot of human time. 

Certain optical methods, such as spectroscopic ellip- 
sometry, or a combination of reflection and transmission 
do not formally rely the KK relations, since both t\ and 



€2 can be obtained directly from the measured data at 
each frequency. Even though the independence on the 
KK relation is a big advantage of such techniques, the 
method, proposed in this paper, can yet be very useful 
here as a stringent test on the overall KK consistency of 
the experimental data. 

In this article we assumed that the materials under 
study are non-magnetic, which means that magnetic per- 
meability fi is equal to 1. For magnetic materials, the 
method can be used to determine n(u>), which also satis- 
fies the KK relations. Of course, in order to disentangle 
n(u>) and eioS) one may need to analyze several measure- 
ments simultaneously (for instance, the transmission of 
samples of different thicknesses) . 

So far we discussed the linear-response optical mea- 
surements only, but, of course, the range of subjects, 
where the KK relations are used, is not restricted to the 
linear optics. The method is potentially applicable to 
any spectroscopic technique, where the measured spec- 
tra depend on a combination of real and imaginary parts 
of a causal, a therefore KK-compliant, response function 
to be determined. The examples include photoemission, 
neutron scattering 45 , X-ray spectroscopy^^, electron 
energy loss spectroscopy (EELS)2I, magneto optical Kerr 
effect (MOKE) 4 «, Josephson-contact spectroscop y 4 ^ 50 , 
electrode impedance measurements^, electron spin 
resonance (ESR)£2iS j low-energy electron diffraction 
(LEED) 54 , acoustics 55 ! 56 and others. Since KK-type rela- 
tions exist for higher-order susceptibilities 57 58 , this way 
to analyze data can also be useful in non-linear spectro- 
scopic techniques. 

The author is grateful to D. van der Marel, F. P. Mena, 
H. J. A. Molegraaf and F. Carbone for numerous discus- 
sions and invaluable practical help. Useful advices about 
solutions of ill-posed inverse problems were given by S. 
V. Rotin. 
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